Mon. Not. R. Astron. Soc. 000,[TH5](2008) Printed 8 October 2008 (MN KTeX style file v2.2) 



Galactic dynamo action in presence of stochastic alpha and shear 

Sharanya Sur\ Kandaswamy Subramanian^ * 

^Inter-University Centre for Astronomy and Astrophysics, Post Bag 4, Ganeshkhind, Pune411 007, India 



ABSTRACT 

Using a one-dimensional aw-dynamo model appropriate to galaxies, we study the possibility 
of dynamo action driven by a stochastic alpha effect and shear. To determine the field evolu- 
tion, one needs to examine a large number of different realizations of the stochastic component 
of a. The net growth or decay of the field depends not only on the dynamo parameters but also 
on the particular realization, the correlation time of the stochastic a compared to turbulent dif- 
fusion timescale and the time over which the system is evolved. For dynamos where both a 
coherent and fluctuating a are present, the stochasticity of a can help alleviate catastrophic 
dynamo quenching, even in the absence of helicity fluxes. One can obtain final field strengths 
up to a fraction 0.01 of the equipartition field B^q for dynamo numbers \D\ ^ 40, while 
fields comparable to B^q require much larger degree of a fluctuations or shear This type of 
dynamo may be particularly useful for amplifying fields in the central regions of disk galaxies. 
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1 INTRODUCTION 

Large-scale magnetic fields in stars and galaxies are thought to 
be generated and maintained by a me an-field turbulent dynamo 
jMoffatll 1 19781 : iKrause & Riidiej |l98d). The potential driver of 
such mean-field dynamos is the a-effect, arising whenever one 
has rotation and stratification in a turbulent flow. Mean-field dy- 
namo (MFD) models using a coherent a-effect and shear have 
been invoked to explain large-scale fields observed in disk galaxies 
( lRuzmaikinetaI.lll988h . 

The possibility of efficient dynamo action arising from ran- 
dom fluctuatio ns in the a-effect in combination with shear was first 
pointed out bv lVishniac & Brandenburj jl997h . They investigated 
a reduced mean-field dynamo model appropriate to accretion disks 
and showed that growth can occur for large enough random fluctua- 
tions in alpha. Several authors have since elaborated various aspects 
of this stochastic alpha-shear dynamos iSokolofl 19971; Silantev 



20001 : iFedotov et a0 2006^: 'Proctoi^2007l: iKleeorin & Rogachevskii 



20081) . In particular, Sokoloft ( 1997) examined a model of a disk 



dynamo with a fluctuating alpha antisymmetric in space but which 
changes sign randomly with equal probability. He argued that inter- 
mittent large-scale magnetic fields can grow. The role of a stochas- 
tic a has also been analyzed in the context of solar dynamos 
JProctor 2007; Brandenburg & Spiegel 2008; Moss et al. 2008). 

The exact origin of such an incoherent a-effect is as yet un- 
clear. In any large Reynolds number system, many degrees of free- 
dom exist, and hence there could always be a stochastic compo- 
nent of the mean turbulent electromotive force (emf). This could 
lead to additive or a multiplicative noise in the MFD equations. 
Additive noise provides a seed field for the dynamo, whereas mul- 
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tiplicative noise in say the a effect, combined with shear, c an lead 
to exp onential growth of the mean field. In the solar context. lHovnd 
( Il993h argued for a fluctuations ~ uo / VM, where uo is the tur- 
bulent velocity and M is the number of cells being averaged over 
in defining the mean field. In principle this can be larger than 
any coherent a-effect. Multiplicative noise is also seen in simu- 
l ations which m easure the a-effect both in t he kinematic regime 
(Suret al. 2008) and in the no nlinear regime ( Cattan eo & HughesI 
2006.; Brandenburg etal .120081) and also in direct simulations of the 

ijalactic dvnamo jGressel et al.ll2008l) . In fact [Brandenburg et al.l 
20081) measure an incoherent a-effect, with a Gaussian probabil- 
ity density function (PDF), even in turbulence driven with a non- 
helical forcing, where one does not expect a coherent a-effect. 
Com bined with shear, such systems show large-s cale dynamo ac- 
tion terandenburg et alJbOOsI: lYousef et alj|2008h . Here, we sim- 
ply examine, in the context of galactic dynamos, the consequence 
of having an incoherent alpha effect, without considering in detail 
its exact origin. The growth of the mean field varies significantly 
from one real ization of the s tochastic process to another, as also 
pointed out in lSokolo3 ( Il997h . It is therefore necessary to examine 
a large number of realizations of the stochastic a to determine the 
efficiency of the stochastic aoj-galactic dynamo. 



We outline in section 2, the basics of a one-dimensional 
stochastic aoj-dynamo model appropriate to galaxies. We present 
numerical solutions of the above model in section 3, with two dif- 
ferent PDF's fo r the stochastic alpha; the first as considered in 
ISokolofd l fT997h and the second where the stochastic alpha has a 
gaussian PDF. In Section 4, we explore the possibility of alleviating 
catastrophic a-quenching in absence of helicity fluxes by including 
the effects of a stochastic a. Section 5 summarizes our results and 
the implications of a stochastic atj-dynamo for galaxies. 
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2 THE STOCHASTIC ALPHA-SHEAR DYNAMO 

In MFD theory, one starts by splitting tlie relevant pliysical quanti- 
ties into mean and fluctuating parts, for example B = B + b for 
the magnetic field and [7 = [/ + u for the velocity field. The over- 
bars denote a suitable averaging procedure with b = u = Q. This 
results in the standard mean field dynamo equation 



dB 

dt 



^ =\/ X {U X B + £ -r]\7 X B), V ■ B = 0. 



(1) 



The averaged equation now has a new term, the mean electromotive 
force (emf) £ — u x b, which crucially depends on the statistical 
properties of the small-scale velocity and magnetic fields, u and b, 
respectively. U is the mean fluid velocity. Assuming that B is spa- 
tially smooth, £ can be expressed in terms of B and its derivative. 



u X b — aB — rjtJ 



(2) 



Here J — V x B/fio (we assume fig — 1 hereafter) and a and rjt 
are turbulent transport coefficients that can be expressed in terms of 
the statistical properties of the flow. In the kinematic regime, and 
assuming isotropic turbulence, one has a — a]^ — — |t?x ■ V x u, 
and the turbulent diffusion coefficient rit = ^tu'^. Here, r, the 
correlation time of the turbulent velocity u, is assumed to be short. 

Since galactic disks are thin, it often suffices to consider a one 
dimensional model, where only z de rivatives of physical variables 
are retained jRuzmaikin et alj|l988l) . For the stochastic dynamo 
which we examine here, we also modify the a-effect to be of the 
form: a = Ok = + ai{z,t). Here a{z) is the average qj;, 
while ai{z, t) is the stochastic a term. Therefore, the total a is a 
sum of the standard kinetic alpha a and a stochastic component ai . 
Further, we consider a mean flow consisting of only a differential 
rotation such that U — (0, rfi(r), 0). Then, going to dimensionless 
variables, Eqn. l[T} gives evoluti on equations for the azimii t hal (B ^) 
and radial (Br) fields, (see also lVishniac & Brandenburg ( Il997h ) 



^ = {Rc.g{^)B, + Q^f{z)N B,) + ^ 



(3) 



dB^ 
dt 



= R^Br 



d_ 

dz 



Rc,g{z)Br + Q, 



'n-)NB,.) + ^.(4) 



Here the length and time units are h and td = /rjt respectively, 
with h the semi-thickness of the disk. We adopt a = aog{z), and 
ai = asf{z)N{t) where f{z) = g{z) — sm{TTz) takes care of 
the symmetry condition. N{t) is a stochastic function. In our nu- 
merical solutions we adopt the following procedure: We split t into 
equally spaced intervals [utc, (n + 1)tc], where Tc is the correla- 
tion time of the stochastic alpha, and n = 0, 1, 2.... are integers. 
And in any such time interval A*' is a random number chosen from 
a Gaussian (or some other) probability distribution, with unit vari- 
ance. The relevant dynamo control parameters are Ra, Qa and Rt^ 
defined as 



_ aph _ ash 



Gh^ 



(5) 



Here G — rdQ,/dr = —Q,, for a flat rotation curve. From Krause's 
formula, qo — I'^O./h, where Zq is the integral scale of interstel- 
lar turbulence. Then Ra ~ Silted, assuming rjt ~ IqUo/?, and 
T ~ ted = /uo, the eddy turnover time. Typical values of the dy- 
namo control parameters in the solar neighborhood of the Milky 
Way could be Ra ~ 1.0, and \R^\ 10 — 15, correspond- 
ing to a " dynamo number" D — RaRui ~ — 10 to D ~ —15 
jRuzmaik in et al. 1988). The strength of Qg ~ 3{h/lo)M-'^^^ if 
one uses the estimate of lHovnd jl993l) . A horizontal average over 
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Figure 1. Time evolution of the rms mean magnetic field in the Sokoloff 
1997 model for different realizations using a long correlation time for 
oi (t). Parameter values used are Ra = 0.0, Qa = 1.0, and iJ^ = —40. 
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Figure 2. Time evolution of the rms mean magnetic field in different real- 
izations using a short congelation time for at (t) and a gaussian PDF. Plots 
are obtained for parameter values Ra = 0.0, Qa = 1-0, and R^ = —40. 



a scale h to define B (cf. lSrandenburg et"aD (l2008l) ) would suggest 
M ~ (h/lo)^ and hence Qa ~ 3. However since the exact ori- 
gin of such fluctuations is as yet unclear, we will vary Qa around 
these values. Thus in general, we will have \Ruj\ 3> -Ra, Qa so that 
one can make the standard atj-dynamo approximation, where one 
neglects the terms with co-efficients Ra and Qa in Eqn. l|4j. 

Note that Q oc 1/r, and thus one can have larger dynamo pa- 
rameters towards the disk centre, depending also on how h and lo 
behave there. The disk height could be smaller, but lo could also be 
smaller in the denser inner galactic regions, where supemovae are 
more confined. This could lead to a net increase in R^j oc 
Any increase in Ra depends on how much /q decreases compared 
to the increase in Q.. Changes in Qa depend on the origin of the 
a fluctuations. For example, if h decreases by factor 2 and /q de- 
creases a factor 5 in the inner galaxy, Rt^ would increase by a fac- 
tor 6.25(r/2kpc)~^ and Ra or Qa would remain about the same, 
compared to the solar neighborhood. Overall larger dynamo num- 
bers can be expected in the inner regions of disk galaxies. We now 
turn to the solution of the stochastic au- dynamo equations. 
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Figure 3. Frequency distribution of tlie dynamo amplification A = B/Bq 
at t = 2 (blue/thin), t = 10 (red/normal) and also at i = 15 (black/thick) 
(in panel d), for 1000 realizations of ai (t) in each histogram. Here Ra = 
0.0, Qa = 1.0, and = -40, -80, -120 and -180. Bin size is 0.2. 

3 NUMEMCAL SOLUTIONS 

Our primary interest is in a scenario where large-scale dynamo ac- 
tion is possible in presence of stochastic alpha and shear. Thus we 
first seek numerical solutions to Eqs. l[3) and l[4j, in the aoj-dynamo 
approximation, with the coherent part of the a-effect taken to be 
zero; that is with Ra — 0. The code uses a 6th order explicit finite 
difference scheme for the space-derivatives and 3rd order accurate 
time-stepping scheme; see iBrandenbur j ( l2003h for details. We use 
vacuum boundary conditions for the fields 

Br = B^ = atz = ±h (6) 

As a test case, we numerically implemented the ISokoloflj 

( Il997h model with a = and ai = asf{z)N{t); N being ei- 
ther + 1 or —1 with equal probability in any time interval utc < 
t < {n + 1)tc. For = 1, the system behaves as a standard 
QLj-dynamo with growing solutions, while for A'^ — —1 we have 
decaying oscillations. So if the system is evolved over a finite time 
interval, there would be random instances of growth and decay. If 
7 is the growth rate of the growing solutions and — that of the de- 
caying ones, the ensemble averaged growth rate is F = (7 — (") /2 
jSokolofdl 19971) . Thus when 7 > C> one obtains F > resulting in 
an overall growth above a critical dynamo numbe r Dc - To estimat e 
Dc, we use the perturbation solutions discussed in ISuretalJ l l2007h . 
This gives 7 ~ — 7r'^/4 + \/ ■ k\D\/2 and <^ ~ 7r^/4, and thus 
for the lSokolofd jl997i) model. This is somewhat larger 
in magnitude than the critical dynamo number ~ 7r^/4, which ob- 
tains for the coherent au dynamo (by demanding 7 > 0). 

These features are illustrated in Fig. (T). Here we have chosen 
Tc = 2td so that one can clearly see both the growing and decay- 
ing phases and their net effect. Starting with random seed fields 
Br, B^ ~ 10~®, we find a number of growing as well as decay- 
ing realizations for moderate dynamo number D = —40. It is ev- 
ident from Fig. (TJ, that any given realization has periods, N+Tc, 
of steady growth (when A^ = 1) and periods, N-Tc of oscillatory 



decay (when A = — 1). One gets a net growth of the field in about 
65% of the realizations, as roughly expected from the above ar- 
guments for \D\ > I -Del- For a larger magnitude of the dynamo 
number \D\ one gets a greater probability for growth. We have also 
examined the opposite limit when Tc < td, and find that the dy- 
namo becomes less efficient (see below). These solutions clearly 
demonstrate the basic idea b ehind the incoherent au dynamo as 
discussed in ISokolofj ( IT997h . that one needs to consider many re- 
alizations of the stochastic process. Just solving a double averaged 
version of the MFD equations need not be representative of the ac- 
tual evolution of the dynamo for a given realization. 

Of course th e PDF of the sto chastic alpha is not expected to 
be as described in ISokolofdfl997h ; for example [Brandenburg et akl 
( I2OO8I) found it can be approximated as a Gaussian. Also in gen- 
eral we expect Tc < td- We now present the results obtained by 
solving Eqs. ([3} and ([4j with the random number A for qi cho- 
sen from a gaussian PDF, and adopting Tc = G.Q2td. This is about 
l.bted, assuming h = 500 pc, rit = 10^^ cm^ s~^, Iq = 100 pc, 
and uo ~ 10 km s~^. The initial seed fields are random with ampli- 
tudes of 0(1). The MFD equations were evolved up to 10 turbulent 
diffusion time scales, t — 10, for dynamo numbers D = —40, —80 
and —120 and upto t = 15 for D = —180. We also considered 
1000 realizations of ai{t) for each D so as to obtain good statis- 
tics. Note that to probe the PDF of the dynamo amplification upto 
a 3(T level one needs about these many realizations. 

Fig. ^ shows the time evolution of the RMS (large scale) 
magnetic field, B, for a subset of realizations with Ra ~ 
0.0, Qa ~ 1.0 and R^j = —40. There is an initial decay of B, 
while the system discovers the proper eigenfunction. Further evo- 
lution then occured on the diffusion time-scale td- In all realiza- 
tions, B shows an oscillatory decay, even though a significant num- 
ber of realizations showed growth of _B up to t = 2. For higher 
dynamo numbers, growth is sustained for a lon ger time and for 
a larger number of realizations. We find that the 'So kolofj h997h 
model also shows similar features for short correlation times. Thus 
having Tc <^ td, qualitatively changes the behavior of the dynamo 
and leaves an imprint of td in the system evolution rather than Tc- 
In order to have a quantitative measure of how many realizations 
show net growth, we show in Fig. [3] the frequency distribution of 
the dynamo amplification A — B/Bq at t = 2, t = 10 and also at 
f = 15 in panel (d) for 1000 realizations of ai (t), at dynamo num- 
bers, D = -40, -80, -120 and -180. Here Bo = 0.32 - 0.35 is 
roughly the value to which B initially decays in all the realizations. 
For D — —40, —80 and —120, we obtain A > 1, for respectively 
34%, 65% and 82.8% of realizations at t = 2. However at a later 
time t — 10, this percentage decreases to 0%, 18% and 24% re- 
spectively. This is evident in the gradual shift of the histogram to 
the left. For \D\ = 160 - 180 the PDF of \A\ remains stationary 
at late times; see panel (d). Above this range, the mean amplifi- 
cation secularly increases with time. Thus, our results show that a 
stochastic aoi-dynamo is reasonably efficient over a few td even at 
Rui = —40, but requires much larger dynamo numbers, as plausi- 
ble towards galactic centres, to sustain fields for long periods. 



4 DYNAMICAL ALPHA QUENCHING OF THE 
STOCHASTIC DYNAMO 

Conservation of magnetic helicity is regarded as a key con- 
strai nt in the evolution of large-scale magnetic fields (see 
Bran denburg & S ubramanian (2005) for a review). A consequence 
of helicity conservation is the production of equal and opposite 
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amounts of magnetic helicity in B and b by tlie turbulent emf £. 
Closure models then imply a suppression of dyna mo action due 
to the growing current helicity as s ociated with b |Pouquet et alj 
Il97d iKleeorin & Ruzmaikin|[l982l : iBlackman & Fieldl 20021) . The 
effect of the small-scale magnetic field on the total a-effect is 
described by the addition of a magnetic alpha to the kinetic al- 
pha, Q = ajj + am- Here a-^ represents the kinetic a-effect 
and Om = Ip'^^tj ■ b is the magnetic contribution to the a- 
effect, with p the fluid density. Specifically, the growth of the 
magnetic alpha (am) to cancel the kinetic alpha (ak) results 
in a suppression of the total a-effect. This suppression can be 
catastrophic in the sense that the large-scale field is quenched 
in an Rm dependent manner. Helicity fluxes across the bound- 
aries of the disk have been identified as a possible mechanism to 
shed small-scale magnet i c helicity, and prevent such quenching 



telackman & Field 2000; Kleeorin et al. 2000: Blackma n & Field 
2001 - Vishniac & Cho, 2001. : .Brandenbura 2QQ5: . .Shukurov et al, 
2006l : ISuretal.ll2007h . 



This situation could change in the presence of a stochastic 
component, as the kinetic alpha can undergo frequent sign rever- 
sals. Hence by the time the am grows to cancel ak, the kinetic 
alpha itself might have changed sign. It is then of interest to ask 
whether addition of a stochastic component to the kinetic alpha can 
stem the catastrophic quenching. This would then naturally pro- 
vide a mechanism for healthy dynamo action even in the absence 
of helicity fluxes. The numerical analysis of the previous section is 
therefore extended by including an am contribution to a in equa- 
tions ([3} and ^ supplemented with an evolution equation for am. 
This can readily be motivated by considering the helicity conser- 
vation equation written in terms of the helicity density x of the 
small-scale magnetic field jSubramanian & Brandenburgl2006l) . 



^ + V ■ F = -2f ■ B - 2r,j ■ 6. 
at 



(7) 



Here F is the helicity flux density. Retaining only the z-derivatives 
and using the fact that the m ain contribution to am comes from the 
integral scale of turbulence dShukurov et~ai]|2006l) , so that 7 ■ b ~ 



fcga • b 



1 ,2 X 

3 p 



(8) 



Eqn. Q can be expressed in dimensionless form, by measuring a 
in units of ao and the magnetic field in units of Beq, where B^^ = 
pu^. In the absence of helicity fluxes, i.e = we have. 



9am 
dt 



-C 



where Rm = Vt/v^ C = 27r^(fco/fci)^, k\ — iv/h and we take 
ki/ko — 5. Further, J ■ B is the current helicity density of the 
large-scale field and is given by 



J ■ rS 04,— Dr—^ — 

OZ oz 



(10) 



We adopt am = at t = and random initial fields of 0(10"^^). 
The system of equations l|4j and l|9) are then solved numerically 
in the atj-dynamo approximation. Note that there is an extra term 
—d{Rc,ainB^) I dz in Eq.|3]and no helicity fluxes are added to the 
r.h.s of Eqn. l|9j. 

Fig. © shows the time evolution of the RMS large-scale 
field in a number of realizations with R^ — 1.0, Qa ~ 0.0 — 
4.0, \R^\ = 40 - 50 and i?m = 10^ Note that for Q„ = 0.0, 
Ra ~ 1.0 and Ru, = —40 (shown by dashed lines), we recover the 
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Figure 4. Time evolution of the rms mean magnetic field for different 
realizations in the dynamical a-quenching model with parameter values 
Re = 1.0, Qc = 0.0 - 4.0, \Ruj\ = 40-50 and = 10^. 
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Figure 5. Space-time diagrams of the radial and azimuthal components 
of the large-scale tield for a realization with parameter values Ra = 
1.0, Qa = 4.0 and \Ru;\ = 50. The color bars on the left panel shows 
the magnitude of the field. 



standard result that the magnetic field is catastrophically quenched 
to very low values. The catastrophic quenching still obtains in some 
realizations for a moderate value of Qa = 1.0 (shown in dotted line 
in the figure). But in other realizations (shown by dash-dotted lines 
in the figure), a stochastic kinetic alpha alleviates this quenching to 
obtain fields of order 0.01 — O.OOlBeq. A detailed analysis shows 
that, for these dynamo parameters, B has a net growth in about 
13% of all the realizations, even till t = 20. In fact, stronger val- 
ues of Qa and Ru, can even amplify the field to near equipartition 
values. Such an example, adopting Qa = 4.0 and \Rui\ ~ 50 is 
shown by the solid line in the above figure. A space-time diagram 
for this realization, between times f = 8 — 14, is shown in Fig.|5] 
Both the radial and azimuthal fields have quadrupolar symmetry 
and show several reversals in sign during this period. We recall that 
high values of the dynamo control parameters are plausible towards 
the central regions of a galaxy. Therefore a stochastic aoj-dynamo 
is more likely to grow coherent magnetic fields efficiently towards 
the central regions of disk galaxies. 
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5 CONCLUSIONS 

We have examined here how a stochastic a-effect in association 
with shear can lead to the generation of large-scale galactic mag- 
netic fields. To determine the field evolution, one needs to examine 
a large number of different realizations of the stochastic Qi (t) . The 
net growth or decay of the field depends on the particular realiza- 
tion, the correlation time of the stochastic a compared to turbulent 
diffusion timescale and the time over which the system is evolved. 

The result s are illustrated first with the simple model of 
ISokolofj j 19971) in Fig. [T] Here the magnitude of qi, takes ran- 
domly a value +a3g{z) or —asg{z) over any time interval Tc, with 
equal probability. Any given realization of ai (t) will have A'^+ pe- 
riods of steady growth (when A'^ = 1) and A'^_ periods of oscil- 
latory decay (when A'^ — —1). But, since the growth rate (7) and 
decay rates (Q ar e differen t , this could lead to a net growth or decay 
as pointed out bv lSokolof3 jl997i) . The critical dynamo number for 
getting growth for say 50% of realizations is \Dc\ ~ 30, moder- 
ately larger than \Dc\ ~ 10 required for the coherent atj-dynamo. 
Our numerical solutions confirm the applicability of this picture 
for long correlation times Tc = 2td, while for Tc < td, the picture 
changes qualitatively, and the dynamo becomes less efficient. 

We then examined more realistic MFD models with a short 
correlation time (tc — LSfed), for a stochastic Qi(t) chosen from 
a gaussian PDF. Our results are given in Fig. [2] and Fig. [3] In this 
case as well, and for a stochastic atj-dynamo number D = —40, 
about 34% of realizations showed growth of B till f ~ 2. However 
subsequently B decays to negligible values by t ~ 10. For higher 
dynamo numbers, growth is sustained for a longer time and for a 
larger number of realizations (cf. Fig.[3]l. One requires D ~ —120 
to obtain long term growth (till t ~ 10), in a significant number 
(« 24%) of realizations, \D\ ^ 160 - 180 for the PDF of \A\ to 
remain stationary and a larger \ D\ for secular growth at late times. 
Note that having an additional coherent a, with \RaRui\ > 10, 
would ensure growth of B. However, the quenching imposed by 
helicity conservation still needs to be alleviated. 

In the usual aoj-dynamo, such helicity conservation leads to a 
growth of the magnetic am, which tends to cancel the kinetic ak, 
so as to catastrophically quench the dynamo. This problem can be 
alleviated by having fluxes of magnetic helicity. In contrast, for a 
stochastic atj-dynamo, by the time Om has grown, Qk could have 
changed sign. This raises the possibility of alleviating catastrophic 
Q- quenching without helicity fluxes. To examine this possibility, 
we solved the stochastic aoj-dynamo equations along with the dy- 
namical a-quenching equation. We included both a coherent and 
incoherent a-effect. In general the radial and azimuthal fields are 
again quadrupolar and show occasional reversals in time (see Fig. 
Is}. When the coherent and stochastic components of a are com- 
parable, and even without any helicity flux, we find that steady 
large-scale magnetic field of strengths of about O.OlBeg could be 
obtained in some realizations even for D ~ —40. Field strengths 
above 0.35^9 are obtained for stronger amplitude of random a- 
fluctuations in association with strong shear (Fig. |4}. Therefore, 
a stochastic acLi-dynamo model is more likely to grow large-scale 
magnetic fields efficiently towards the central regions of a galaxy 
and even in the absence of helicity fluxes. 

We have focussed on the application of a stochastic auj- 
dynamo model to galaxies. However, our emphasis on examin- 
ing a large number of realizations and our results on alleviating 
a-quenching with a random a could be applicable to other astro- 
physical dynamos as well. More work is needed to elucidate the 



origin of the incoherent a and also study the influence of spatial 
decorrelation on the dynamo. 
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